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ABSTRACT 

We present the first high-angular resolution study of the MonR2 star-forming com¬ 
plex carried out with the Submillimeter Array at (sub-)millimeter wavelengths. We 
image the continuum and molecular line emission toward the young stellar objects in 
MonR2 at 0.85 mm and 1.3 mm, with resolutions ranging from 0.5” to ~ 3”. While free- 
free emission dominates the IRS1 and IRS2 continuum, dust thermal emission prevails 
for IRS3 and IRS5, giving envelope masses of ~0.1-0.3 M 0 . IRS5 splits into at least 
two sub-arcsecond scale sources, IRS5B and the more massive IRS5A. Our 12 CO(2-l) 
images reveal 11 previously unknown molecular outflows in the MonR2 clump. Com¬ 
paring these outflows with known IR sources in the IRS5 and IRS3 subclusters allows 
for tentative identification of driving stars. Line images of molecular species such as 
CH 3 CN or CH 3 OH show that, besides IRS3 (a well-known hot molecular core), IRS5 is 
also a chemically active source in the region. The gas excitation temperature derived 
from CH 3 CN lines toward IRS5 is 144 ± 15 K, indicating a deeply embedded protostar 
at the hot-core evolutionary stage. SED fitting of IRS5 gives a mass of ~7M 0 and a 
luminosity of 300 L 0 for the central source. The derived physical properties of the CO 
outflows suggest that they contribute to the turbulent support of the MonR2 complex 
and to the gas velocity dispersion in the clump’s center. The detection of a large num¬ 
ber of CO outflows widespread across the region supports the competitive accretion 
scenario as origin of the MonR2 star cluster. 
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Introduction 


Massive young stars are found almost exclusively in clusters with stars of lower stellar masses. 
How molecular clouds collapse and fragment to assemble stars of distributed masses has been the 
subject of intensive studies in the past several years. There are currently two main theories of mas¬ 
sive star and star cluster formation. While the first scenario involves the monolithic grav itational 


collapse of a turbulent-supported massive core a t high accretion rates (e.g. lYorke fc Sonnhalter 


2002 


McKee fe Tan 1120031 : iKrumholz et al. Il2009f 'l. the second scenario - the competitive accretion 


scenario - initially forms a cluster of low- and intermediate-mass stars where the most m assive stars 


grow via Bondi-Hoyle accretion at the center of the cluster potential well (see e.g. Bonnell fc Bate 


20061 1. In the competitive accretion scenario, it is thus expected that massive stars form at the 
center of the cluster with the stellar density of low-mass stars peaking at this position. At the 
cluster center where stellar densities are predicted to be high, dynamical encounters may take 
place as proposed f or the Orion BN/KL and the D R21 high-mass star forming regions (see e.g. 
Zapata et al. II20091 . 120131 : iRivilla et al. |[2013a| , [201J ). Therefore, to study the formation of massive 
stars, it is crucial to better understand the distribution of low-mass stars and their impact on the 
subsequent evolution of high-mass star forming clusters. 

Molecular outflows appear early in the star formation process as a mechanism to remove 
the excess angular momentum in the star/circumstellar disk system produced by gas accretion. 
Molecular outflows thus represent a useful tool to probe the population of Class 0 low-mass stars 
in clusters. Since molecular outflows inject momentum and energy into the surrounding medium, 
they are expected to increase the level of turbulence and even disrupt the parental m olecular cloud 


potentially impacting any subsequent star formation within the cluster (see IStanke fc Williams 
2007; Rivilla et al. 2013b). 


The Monocero s R2 (MonR2) molecular cloud complex is located at a distance of 830 pc 
(|Herbst &: Racine 1119761 ) and e xtends over 3°x6° in the pl ane of the sky. MonR2 contains sev¬ 
eral sites of star formation (see ICarpenter Hodapp 120081. and referenc es therein). The MonR2 
clump (size of ~3' or 0.7 pc) has a mass of 1800 M 0 (Ridge et al. 2003) and is the m ost active 


star f orming region in the MonR2 complex as revealed by th e blister-type HII region (IMassi et al. 


1985 


Wood fc Churchwell 1 19891 ) , the dense stellar cluster (jCarpenter et al. 119971 : iKohno et al. 


2002m the large -scale CO outflow (one of the largest - ~0.6-7 pc - and most massive - 275 M 0 


reported so far; Loren 


found in this region (Downes et al 


1981 


Wolf et al. 


19 


3 


1990; 


Tafalla et al. 


Genzel fc Downes 


IR sources detected in the cluster (Car pente r et ah 19971 b the sources IRS1, IRS2 and IRS3 are the 


1977; 


199 7b and the HoO and OH masers 


Smits et al. 


1998). Among the 


most luminous objects in the Mo nR2 region, with luminosities of 3000 L 0 , 6000 L 0 and 15000 L 0 
respectively ([Henning et al. Ill992l h The ionizing source of the H II region is found 4” away from 
the peak of the radio continuum emission (jAspin Sz Walther 1119901 ) and it is resolved into two com - 
ponents, IRS1 NE and IRS1 SW. While IRS1 NE is a foreground field star ( Howard et al. 1994 b 


IRS1 SW is responsible for the HII region and it is likely a BO-type star in the ZAMS (jMassi et al 


1 98 5). IRS2 appears as a bright source in the K-band, indicating that this object is embedded 
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in a molecular cloud (jBeckwith et al. I Il976l : lAspin fe Walther II1990T). IRS2 is very compact and 


2013 

). This source is respc 

msible for the reflection nebula seen in the K-band and nbL images 
l coincident with the edges of the radio continuum emission mapped 

(Carpenter et al. 

1997) anc 

with the VLA 

Massi et al. 

1985 

). The IRS3 source is the most luminous object in the region and, 


as revealed by Speckle^interferometry, it is in fact a cluster of IR sources, one of which (IRS3B) 
shows a microjet (Preibisch et al. 2002). Since the orientation of this jet (position angle of 
P.A.~50° from north to east) is roughly perpendicular to the large-scale CO outflow, it is unlikely 
that this source (IRS3B) is the driving object of the outflow. In fact, it still remains unknown 
which so urce within the MonR2 clust e r is the actual drivi ng source of the 7pc-long CO outflow in 
MonR2 ( Giannakopoulou et al. 1997 : Tafalla et al. 1997 1. 


Other less luminous IR sources are also detected toward the MonR2 cluster, for example 


IRS4-7, aN, as, and b-h (see iBeckwith et al. Ill976l : lHackwell et al. Ill982l : lAspin fe Walther lll99Cl : 


Carpenter et al. 1997). Physical properties have only been constrained for a few of them: to- 
tal luminosity estima t es of ^700-800 a nd ~300 L 0 exist for IRS4 and IRS5, respectively 


(Henn ing et al. 


1992 


H ackw ell et al. 1 982 ). More extensive studies, especially at high-angular 


resolution and across several wavelengths, are needed to fully characterize the properties of these 
intermediate-luminosity sources associated with the MonR2 cluster. 


By using near-IR 2.12 /xm H 2 observations, Hodapp (2007) reported the detection of 15 new H 2 
jets in the MonR2 cluster that are likely associated with young Class 0 and I sources at their main ac¬ 
cretion phase. These jets were mainly found toward the edge s of the MonR2 molecular cloud, which 


was interpreted as a signature of triggered star formation (Hodapp 120071 : ICarpenter fc Hodapp 


20081 ). However, near-IR observations may provide a biased view of the global population of molec¬ 


ular outflows in deeply embedded star forming clusters such as MonR2, because of the large ex¬ 
tinction found in these regions. Alternatively, the rotational transitions of carbon monoxide (CO) 
in the millimeter and sub-millimeter wavelength range are known to be excellent probes of the 


Gueth & Guilloteau 

1999; 

Lee et al. 

2007) 


Interferometric observations of the rotational transi¬ 
tions of CO at sub-millimeter/millimeter wavelengths are thus needed to unveil the most embedded 
(and the youngest) population of molecular outflows, and therefore of Class 0 and Class I low-mass 
stars, at the cluster centers of high-mass star forming regions. 

In this paper, we report the first interferometric images of the therma l continuum and molecular 
line emission carried out with the Submillimeter Array (SMA; Ho et al. 2004 *1 toward the MonR2 
high-mass star forming cluster at angular resolutions ranging from 0.5” to ~3” (i.e. from 0.002 to 
0.012 pc at a distance of 830pc)@. Our SMA images of the 12 CO J=2—M line emission reveal the 


1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 

2 MonR2 has only been imaged once before in HCO+ J=l—by using the Hat Creek interferometer and at an 
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presence of 11 new CO outflows toward the innermost ~l , -region in the cluster. This population 
of outflows is widely distributed across the MonR2 clump, which contradicts the proposed scenario 
of triggering for the MonR2 star cluster. The paper is organized as follows. In Section 2, we report 
the details about the SMA observations. The results of the dust continuum and 12 CO line emission 
are presented in Sections 3.1 and 3.2, respectively. In Section 3.3 we analyze the line transitions 
of other molecular tracers such as CH 3 CN, CH 3 OH or SO 2 detected toward the IRS3 and IRS5 
sources, and in Section 4 we discuss our results. The conclusions are summarized in Section [5j 


Observations and data reduction 


Observations of the MonR2 star forming region were carried out with the SMA at 1.3 mm in 
very extended configuration (VEX) in February 2010 for two tracks, and in Compact configuration 
in November 2011 for one track. Because the Compact track was performed in dual receiver 
mode, we simultaneously imaged MonR2 at 0.85 mm. The observations were ph ase-centered o n 


1997 

), and the central radial velocity of the source was set at 10 km s 1 

Torrelles et al. 

1983) 


Instrumental and calibration parameters for the VEX and Compact observations are summarized 
in Table [B The VEX observations were obtained in single-receiver mode with a total bandwidth 
of 4 GHz per sideband. The resulting frequency coverage was 216.64-220.60 GHz and 228.60- 
232.58 GHz. For the Compact observations, the receivers were configured in dual receiver mode 
with a total bandwidth of 2 GHz per sideband. As a result the frequency coverage was in 2 GHz 
windows, as follows: 220.24-222.23 GHz and 230.24-232.23 GHz for the 230 GHz receiver, and 
342.67-344.66 GHz and 352.67-354.66 GHz for the 400 GHz receiver. The high-frequency receiver 
on one antenna was malfunctioning during the Compact track and the resulting data flagged, 
causing the beam to be elongated at 400 GHz. The correlator provided a uniform channel spacing 
of 0.8 MHz, yielding a channel width of 1.1 km s 1 at 1.3 mm (230 GHz) and ~ 0.7 km s 1 at 
0.85 mm (400 GHz). The primary beam (field of view) of the observations was ~57” at 230 GHz 
and ~35” at 400 GHz. Calibration of the raw data was done with the IDL MIR software package, 
while continuum subtraction, imaging and deconvolution were carried out with MIRIAD. We note 


angular resolution of ~9” (|Gonatas et al.lll992f ). 


Table 1. Instrumental parameters of the SMA observations. 


Date 

Config. 

L.O. freq. (GHz) 

Synthesized Beam 

T225GHZ 

Tsys (K) 

BP Cal. 

Flux Cal. 

Gain Cal. 

Feb 13, 2010 

VEX 

224.611 

0.68" x 0.45", P.A.=-43° 

0.05 

200-240 

3C273 

Titan 

0607-085/0730-116 

Feb 20, 2010 

VEX 

224.611 

0.68" x 0.45", P.A.=-43° 

0.03 

200-240 

3C273 

Titan 

0607-085/0730-116 

Nov 15, 2011 

Compact 

226.227 

2.75" x 2.73", P.A.=-77° 

0.07 

200-240 

BLLAC 

Ganymede 

0607-085/0530+135 

Nov 15, 2011 

Compact 

348.324 

2.29" x 1.50", P.A.=-39° 

0.07 

400-500 

BLLAC 

Ganymede 

0607-085/0530+135 
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that for the Compact data, the quasar BLLAC was used as bandpass calibrator (see Tabled]). 
Since BLLAC’s spectrum shows a prominent 12 CO absorption feature, the bandpass calibration of 
the corresponding upper sideband was performed by phase-conjugating the lower sideband SMA 
BLLAC data. 


3. Results 


3.1. Dust Continuum Emission 


3.1.1. Compact SMA Observations 


In Figure dl we present the continuum emission seen in Compact configuration at 1.3 mm 
(left panel) and at 0.85 mm (right panel) towards MonR2. For com parison . Figure dl includes the 
1.3 cm Very Large Array continuum image from the dataset of Zapata et al. ( 2009h . where the 
MonR2 central HII region is clearly visible. Our SMA image at 1.3 mm shows extended structure 
toward IRS1 coincident with the brightest emission seen at centimeter wavelengths. However, in 
the 0.85 mm map of Figure 1 this source is not detected. This could be due to: i) IRS1 is located at 
the edge of the primary beam of the 0.85 mm observations; and ii) the continuum emission detected 
in the SMA images probes the innermost regions around IRS1 and therefore arises from ionized 
gas associated with the UC HII region. Although subject to large uncertainties, the comparison 
between the 1.3 cm and 1.3 mm flux for IRS1 yields a negative spectral index a~—0.5 (with 
S u cx v 01 ), consistent with free-free emission. By using this spectral index, we can extrapolate 
the expected flux for this source to 0.85 mm. This flux, at the half-power of the primary beam, is 
~20mJy beam -1 , well below the 3cr noise level of 30 mJybeam -1 in the 0.85 mm image of Figured] 
This explains why the IRS1 source has not been detected in our SMA continuum image at 0.85 


Table 2. Properties of the dust continuum emission in Compact configuration. 


Source 

Wavelength 

a (J2000) 

S (J2000) 

Size a (”x”) 

P.A. O 

Fpeak h (Jy/beam) 

F int b (Jy) 

IRS2 

1.3 mm 

06 h 07 m 45.803 s 

-06 o 22'53.51" 

3.15 x 2.78 

75 

0.154 ±0.005 

0.180 ±0.005 

IRS2 

0.85 mm 

06 h 07 m 45.800 s 

-06°22 , 53.55" 

3.00 x 2.79 

140 

0.139 ±0.004 

0.154 ±0.005 

IRS5 

1.3 mm 

06 h 07 m 45.588 s 

-06 o 22'39.42" 

3.14 x 2.79 

122 

0.109 ±0.005 

0.127 ±0.006 

IRS5 

0.85 mm 

06 h 07 m 45.612 s 

—06°22 , 39.47" 

2.95 x 2.86 

10 

0.389 ±0.012 

0.436 ±0.013 

IRS3 

1.3 mm 

06 h 07 m 47.811 s 

-06 o 22'56.18" 

3.37 x 2.88 

55 

0.059 ±0.008 

0.077 ±0.011 


a These are non-deconvolved source sizes. 

b The errors in the peak flux and total density flux correspond to those obtained from the two-dimensional Gaussian fits of 
the emission. The 0.85 mm image was spatially smoothed to the same angular resolution as that of the 1.3 mm image. 
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Fig. 1.— Left panel: Continuum emission detected at 1.3 mm in Compact configuration toward the 
MonR2 sta r forming region (contours) overlaid on the 1.3 cm continuum image (from the dataset 
of Zapata et ah 2009). Contours start at 12 mJy beam" 1 (3a) and are spaced at 24 mJy beam" 1 
(6cr) intervals. Dotted contours correspond to the negative 3cr level. The beam size is shown as a 
red ellipse in the lower left corner. For the VLA data, the first contour is at 3 mJy beam" 1 (3<r) 
and the step level is 17 mJy beam" 1 (17a). The dashed circle indicates the primary beam of the 
observations, ~57” at 1.3 mm. Right panel: As for the left panel, but for the 0.85 mm continuum 
image. The first contour is at 30 mJy beam" 1 (3<r) and the step level is 60 mJy beam" 1 (6a). At 
0.85 mm the primary beam of the observations is ~35” (dashed circle). The plus signs indicate 
the positions of the main sources in the region, including IRS2, IRS3, IRS4, IRS5 and the ionizing 
star of the HII region, IRS1 SW (note that this is differenRfrpm IRS 1 NE, which is an unrelated 
non-embedded star detected in the optical; see Cohen & Frogel 19771 ). The maps were made using 
uniform weighting. 
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From Figure 1, we find that the brightest sources in our SMA images are IRS2 and IRS5. A 
qualitative comparison of the two maps at 1.3 mm and 0.85 mm reveals that these two sources are 
very different in nature. While IRS2 appears brighter at 1.3 mm than at 0.85 mm (as for IRS1 and 
as expected for an HII region), IRS5 shows the opposite behavior, indicating that its emission is 
mainly due to dust. In Table [2] we report parameters calculated from two-dimensional Gaussian fits 
to the continuum emission of these two sources. The derived coordinates are close to those obtained 


from previous infrare d and X-ray measurements (e.g. ICarpenter et, al. 119971 : iKohno et al. 


2ooa 


Nakaiima et al. 112003 1. As IRS5 is offset from the phase center of the observations by approximately 
15”, primary-beam corrected maps were used to estimate its flux parameters. After smoothing the 
0.85 mm images to the same angular resolution as that of the 1.3 mm data (~ 3”), we compare the 
peak flux values at 1.3 mm and 0.85 mm (see Table [2]) to estimate the spectral index, a, of the two 
sources (with S u oc u a ). We obtain cqRS 2 ~ —0.2 and qtrrk ~ 3 .0. The decreasing spectral index 
of IRS 2 is consistent with that previously reported by Jimenez - Serra et al. (J2013I), and indicates 
that this source is dominated by optically thin free-free emission. This source, with an es tima ted 


centr al mass of 13 M 0 (a Bl-type star in the ZAMS) an d a luminosity of 6000 L^ fjHenning et al. 


1992l ). likely powers a dense and collimated ionized jet (|.Timenez-Serra et al. 1120131 ) 


In contrast to IRS2, the spectral index of IRS5 steeply increases with frequency (ams 5 = 3.0), 
indicating that its emission at millimeter wavelengths is dominated by dust thermal emission. This 
suggests that IRS5 is at an earlier stage of evolution than IRS2. Since the emission from IRS5 is 
domina ted by dusk we can estimate the envelope dust mass in the source following the method 
from iHildebrand 1 1 983 ): 

F{v)d 2 


A/dust = 


(1) 


By (Tdust) l^v 

where F(v) is the continuum flux density at frequency v, d is the distance to the source, B u (T^ nst ) 


is the Planck function for dust temperature Tdust, and k„ is the dust opacity. k u (cm 2 g -1 ) is given 
by 10 ( i 2 THz ) ^ ’ where (3 is the dust opacity index. Assuming a power law of the type Q(v) oc i/P 
for the emissivity Q(v), in the Rayleigh-Jeans limit (3 and a are then related by f3 = a — 2 . Since 
qtr.sk ~ 3. 0, we_therefore obtain /3irs5 ~ 1. Although this value is similar to those assumed by 


Qiu et ah (2007) for a sample of massive star-forming regions, we cannot rule out the possibility 
that our estimate of (3 for IRS5qs affected by a small contribution from free-free emission, especially 
at 1.3 mm (see Giannakopoulou et ah 19971 1. As an upper limit for the dust temperature, we use 
the gas temperature derived toward IRS5 from the rotational diagram of the CH 3 CN (12-11) K = 
0 to K = 6 lines detected with the SMA (T ~ 150 K; see Section T3.3D . As a lower limit for thejiust 


tempe rature, we assume Tdust = 50 K, consistent with the Tdust values der ived bv lThronson et ah 


(1980) towards IRS5 and with those estimated by Sridharan et ah (120021 ) towards high-mass star 
forming regions such as IRAS18264-1152 and IRAS23151+5912. 

From the integrated fl ux reported in Table [2l and c onsidering a constant gas-to-dust mass 
ratio of 100 (following e.g. I Giannakopoulou et ah I Il997l : iBeuther et ah l2007blf . we obtain esti¬ 


mates for the IRS5 envelope gas mass of ~ 0.08 — 0.3 M 0 , as reported in Table [3l We note that 
these masses are smaller than those reported at poorer angular resolution (16”-27”) for IRS5 by 















































Giannakopoulou et al. (1997, about 5 Mg, see clump C5 in their Table 4). This is due to the 
fact that the SMA only probes the densest regions toward IRS5, filtering out most of the extended 
emission around this source. In Table [2] we also report the continuum flux measured at 1.3 mm 
toward IRS3 (after the primary-beam correction of the images). Since this source falls outside 
the primary beam of our observations, the measured 1.3 mm flux and derived dust continuum 
properties (Table El) should be taken with caution. 

Toward the north-east of IRS5 we also find some continuum emission in the SMA 1.3 mm 
continuum image at the 9<r level (Figured]). Although the bulk of the emission does not coincide 
with the location of any IR source, the elon gated asymmetry towar d the east of this emi s sion is 


indeed associated with source aN detected by ICarpenter et al. (119971 ) and lAndersen et al. (|2006l) 


As shown below, this source is resolved into several condensations (see Figured]) and could therefore 
represent a cluster of deeply embedded protostars. Future observations covering this region will 
help to establish the origin of the 1.3 mm continuum emission toward the northeast of IRS5. 


3.1.2. VEX Observations of the Small-Scale Structure 


Figure [2] presents the 1.3 mm continuum map of the MonR2 cluster obtained by combining 
the SMA data at 1.3 mm in the Compact and VEX configurations. The VEX data provide a 
much improved angular resolution of 0.5”. Toward IRS1, the extended emission associated with 
this infrared source shows a complex structure with two bright emission peaks with peak intensities 
detected at the 9 a level. These p eaks, however, do not co incid e with the IR sources IRS1 NE or 
IRS1 SW previously reported bv lAspin fc Walther (j 199(1 ) and iHoward et al. (119941 ) and shown 
by black cr osses in Figured! rather, they appear associated with the IR ring nebula related to the 
HII region ( Carpenter et al. 1997 : Andersen et al. 2006; de Wit et al. 2009). This indicates that 
the 1.3 mm emission is likely associated with dense inhomogeneities in the ionized gas. 

Figure [2] also shows that IRS2 remains unresolved by the present 0.76” x 0.54” beam. The 
secondary continuum peak toward the southwest of IRS2 is not associated with a real source, but 
rather due to limited UV sampling. This implies that the source does not include substructure at 
scales of 0.5” (i.e. 400 AU at a distance of 830 pc). This compact morphology is consistent with 
the results of Alvarez et ah (2004) obtained by using near-IR Speckle imaging. 


Further north, IRS5 shows hints of substructure not seen previously in the Compact data alone. 
In particular, a secondary peak to the southeast of the main source appears as a 5cr detection. For 
the purpose of characterizing these sources, we refer to the main IRS5 peak as IRS5A, and to 
the secondary source as IRS5B. Flux parameters for these two sources are reported in Table El 
Following the method in Hildebrand (|l_983j) and assuming a dust opacity index of [3 ~ 1, we 
estimate envelope gas masses of 0.07-0.24 M 0 for IRS5A and of 0.01-0.03 M 0 for IRS5B (see Table 
El) As expected, IRS5A is the dominant source, as it accounts for ~ 90% of the total gaseous 
envelope mass in IRS5 measured in Compact configuration. In addition, an elongation along the 
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Table 3. Envelope mass estimates from COM data. 


Source 

Mass (Mq) 


T du st = 50 K 

T dus t = 150 K 

IRS5 

0.26 

0.08 

IRS3 

0.15 

0.05 


Note. — For IRS3, no emission is de¬ 
tected at 0.85 mm, so the mass is derived 
using the same dust opacity index as for 
IRS 5. 


Table 4. IRS5 properties with combined Compact and VEX data at 1.3 mm. 



IRS 5 A 

IRS5B 

R.A. (J2000) 

06 h 07 m 45.60 s 

06 h 07 m 45.65 s 

DEC (J2000) 

-06°22'39.37" 

-06°22'40.07" 

Source size (”x”, PA) 

0.90x0.70, 44° 

F- 

o 

VI 

Apeak (Jy/beam) 

0.080 ± 0.002 

0.015 ±0.003 

Aint (Jy) 

0.124 ±0.003 


Mass (Tdust = 150 K) 

0.07 Mq 

0.01 Mq 

Mass (Tdust = 50 K) 

0.24 M e 

0.03 Mq 


Note. — The errors in the peak flux and total density 
flux correspond to those obtained from the Gaussian fits 
of the emission. For IRS5B the peak intensity was used to 
derive a mass estimate because the source is point-like and 
no reliable Gaussian fit was obtained. The source sizes are 
not beam-deconvolved. 
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Fig. 2.— Combined VEX and Compact data at 1.3 mm. Left panel: In this map of the entire field 
of view, contours start at 3cr and are spaced at 3cr intervals, with a ~ 2.5 mJy/beam. The beam 
size of the map is shown as a blue ellipse in the bottom left corner of the image. Large crosses 
indicate the position of the main sources in the region as in Figure [H For the top two zoom-in 
panels, we show the 3 a, 4<r, 5a, 8a, 15 a, and 25er contours. For IRS2, we show the 8a, 6a, 10a, 
30(7, 50(7, and 70cr contours. Finally, for IRS1 we show contours starting at 3a and spaced by lcr. 
Dotted contours correspond to the negative 3<7 level. The map was made using uniform weighting. 
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Fig. 3.— Maps of the combined VEX and Compact 1.3 mm data for IRS5 produced with different 
beams. The beam size of the map is shown as a red ellipse in the bottom left corner of each panel. 
The elongation along the northeast-southwest axis of IRS5 remains when the image is convolved 
with a circular beam. We show the 3cr, 4a, 5a, 8a, 15a, and 25a contours, with a ~ 2.5 mJy/beam. 
The map was made using uniform weighting. 
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northeast-southwest axis is highlighted by the contours, suggesting the existence of a third source 
that remains unresolved at the current resolution. As shown in Figure [3] this elongation is robust 
when the image is convolved with a circular beam. Finally, the continuum peak to the west of 
IRS5A might be another source. However, since the emission is only detected at the 4 sigma level, 
higher sensitivity observations are required for confirmation (Figure [3] upper panel). 


3.2. 12 CO Emission: Outflowing Gas in MonR2 


3.2.1. Morphology 


Figure H] presents the integrated intensity maps of the 12 CO(2-l) emission (at 230.538 GHz) 
measured toward MonR2 in compact configuration and overlaid on the continuum emission at 
1.3 mm. The 13 CO(2-l) line transition was also observed within our frequency setup and it will be 
used in Section f3.2.3l to determine the optical depth of the 12 CO outflow emission. The velocity 
ranges chosen for the 12 CO(2-l) emission highlight at least six (labeled with letters in Figured]) 
and up to eleven previously unknown outflows in the central region of the MonR2 cluster. The 
orientations of the outflows are indicated by black arrows in Figure [4] IRS3 and IRS5 are both 
associated with bright and highly collimated bipolar outflows (left panel of Figure (H) The CO 
emission around IRS5 is aligned along the northeast to southwest axis, exhibiting a large degree of 
collimation (outflow B). The redshifted component of the outflow appears between 26 and 40 km s’ 1 
and features two separate 12cr and 15cr peaks. These knots, or ’’bullets”, are commonly observed 
in low-mass outflows such as e.g. L1448-mm and HH211 (see 


2006; 

Palau et al. 

SD 

O' 

o 

CM 

Lee et al. 

2007; 

Hirano et al. 


1992; 


Hirano et al. 


201 


Guilloteau et al. 

Q) and may represent outbursts du e 


to disk variability caused by episodic, unsteady accretion (|Qiu fe Zhang 1120091 : iLee et al. 1 120101) . 
Another redshifted emission peak is detected further to the northeast (+13”, +40”) and appears 
more prominently in primary-beam corrected maps, where it is well-aligned with the rest of the 
redshifted lobe. On these grounds we calculate the physical parameters of the outflow in Section 
13.2.31 assuming all three bullets belong to the redshifted lobe. Its blueshifted counterpart appears 
for both velocity ranges between -16 and -5 km s' 1 and -5 and 1.9 km s’ 1 , and it is less extended 
spatially. IRS5 is located at the geometric center of this outflow and is therefore likely its driving 
source (see Section l3.2.2l) . Similarly, IRS3 is associated with the geometric center of a bipolar 
outflow along the northeast to southwest axis (outflow A). The red- and blueshifted components of 
this outflow are clearly seen in both panels of Figure |4j The blue lobe features two separate knots 
and is spatially more extended. 


The right panel of Figure |4] presents additional outflows identified in the 21 to 26 km s’ 1 and 
-5 to 1.9 km s' 1 velocity ranges for redshifted and blueshifted emission, respectively. A prominent 
large scale redshifted lobe is present along the northwest-southeast direction in the vicinity of 
IRS5 (outflow D). This outflow lacks an obvious blueshifted counterpart. We propose that the 
most likely scenario is that the corresponding blue lobe could be located outside the field of view 
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Fig. 4.— Moment-zero maps of the CO(2-1) transition in MonR2 (colored contours) overlaid on the 
1.3 mm continuum (greyscale). Left panel: The redshifted emission has been integrated over the 
26-40 km s -1 velocity range (a ~ 0.15 Jy/beam km s~ x ); while the blueshifted emission has been 
integrated over the velocity range -16 to -5 km s _1 (a ~ 0.14 Jy/beam km s^ 1 ). The first contour 
and step levels are 4a. Right panel: The redshifted gas appears for velocities comprised between 
21 and 26 km s -1 (a rsj 0.60 Jy/beam km s x ). The 3cr, 5<r and 7a contours are shown; contour 
spacing is 3a thereafter. The blueshifted counterparts appear between -5 km s^ 1 and 1.9 km s” 1 
(a ~ 0.41 Jy/beam km s _1 ). Here the 3cr, 5a and 10 a contours are shown; contour spacing is 5a 
thereafter. The brightest outflows are labeled with letters. Magenta triangles indicate the positions 

(2 0061), while black squares mark the 
(1997). Black arrows indicate the 


Andersen et al. 


Carpenter et al. 


of the IR sources observed with NICMOS by 
coordinates of the IR sources identified by 
conjectured orientation of the outflows. The beam size is shown as a red ellipse in the lower left 
corner, and the dashed circle indicates the primary beam of the observations, ~57” at 1.3 mm. 
The maps were made using uniform weighting. 
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of our observations. The extended redshifted lobe features three approximately equally-spaced 
peaks. The angular separation between the knots is roughly 20”, giving a projected separation of 
approximately 1.5 x 10 4 AU ~ 0.08 pc. If these bullets are associated with periodic outbursts, this 
separation between knots implies time-scales of ~5000 yrs between bursts, which are con sis tent 


with those predicted b y episodic accretion models in low-mass protostars (see iBaraffe &; Chabrier 


201C; 


Zhu et al. 2010). We also detect one red and two blue lobes with aligned elongation axes 


on either side of IRS1 (outflow E). Between IRS3 and IRS5, the image reveals the presence of two 
bright blueshifted condensations (C and F). The emission is strong for a wide range of velocities, 
lacks a well-defined directional axis and well-defined redshifted counterparts. Outflow C is also 
seen in the left panel of Figure|4] for the velocity range from -16 to -5 kms -1 . Finally, we detect 
a number of relatively bright structures located beyond the edge of the primary beam, which we 
tentatively assign to outflows according to their association with sources detected in the near-IR 
across the cluster fSection l3.2.2D . These outflows can be found toward the northern edge of the 
image, west of IRS5, and due north of IRS3 (right panel of Figured). They are unlikely to be density 
enhancements of the large-scale outflow observed by IWolf et al, (1990! ) because there are blue- and 
redshifted lobes present in both the eastern and western parts of the region, while the lar ge C O 


outflow is blueshifted towards the northwest and redshifted in the southeast (ITorrelles et. ah 1983) 


3.2.2. Comparison with IR sources 


The inner region of the MonR2 stellar cluster, approximately the field of view of our observa¬ 
tions, is known to harbor a centrally peaked gas density distribution where the H 2 volum e den sity. 


n(Ho ), peaks toward the position of the IRS 1 and IRS2 sources (see Figures 8 and 9 inlChoi et al. 


the X-rays ([Kohno et al. 


2002 


2000). In a ddition, surveys in the infrared ([C arpent er et al. Ill997l: lAndersen et al. 120061 ) and in 


Nakajima et al. 


2003 


Preibisch et al. 


2002) have also revealed a 


dense population of embedded low-mass stars. Here we compare the clus t er of molecular outflows 


descri bed i n Section 13.2.11 w i th th e IR sources found by ICarpenter et al. (|l997l ) , IPreibisch et al. 


(2002) and Andersen et al. (20061), in order to identify the possible driving sources that power the 
CO outflows detected in MonR2. 

In Figure[5j we present zoom-in maps of the CO outflows detected toward IRS5 and IRS3 (out¬ 
flows B and A in Figure 01) From the left panel of Figure [H we find that the collimated CO outflow 


toward IRS5 is likely associated with the star detected in the near-IR only by Carpenter et al. 


(1997) and which falls close to the location of the IRS5A source identified in the VEX data, the 
most massive star in the IRS5 sub cluster (its envelope mass is the largest one measured in the 
cluster; see Figure [2] and Table 01). The candidate source is marked by an arrow in Figure [5j 

The right panel of Figure [5] shows the pres ence of a stellar subcluster in IRS3 (sources A-F) as 


detec ted from IR data (jPreibisch et al. 1120021 ). As revealed by K-band images (IPreibisch et al 


2002), the massive star IRS3 B presents an elongated nebulous feature that points toward the 
north-east with a P.A. of ~50°. This feature could be the near-IR counterpart of the blue-shifted 
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Fig. 5.— Moment-zero maps of the CO(2-1) transition in MonR2 (colored contours) overlaid on the 
1.3 mm continuum (greyscale). Left panel: A zoom-in on outflow B near IRS5 where the redshifted 
emission has been integrated over the 26-40 km s^ 1 velocity range (a ~ 0.15 Jy/beam km s” 1 ), 
and the blueshifted emission over the -16 to -5 km s^ 1 range (a ~ 0.14 Jy/beam km s^ 1 ). The 
first contour and step levels are 4a. The location of the candidate driving source is marked by a 
black arrow. Right panel: A zoom-in on outflow A near IRS3 where the redshifted gas appears for 
velocities comprised between 21 and 26 km s^ 1 (a ~ 0.60 Jy/beam km s -1 ). The 3<r, 5 a and 7 a 
contours are shown; contour spacing is 3 a thereafter. The blueshifted counterparts appear between 
-5 km s _1 and 1.9 km s^ 1 (a ~ 0.41 Jy/beam km s^ 1 ). Here the 3a, 5 a and 10a contours are 
shown; contour spacing is 5 <t there afte r. M ag enta triangles indicate the positions of the IR sources 
o bserved with NICMOS b y lAndersen et al. (|2006l ) , black squares the stellar pop ulation described 


m 


Car penter et al. (19971) , and black dots and letters the IR sources identified bv lPreibisch et al. 
(2002). The dashed circle indicates the primary beam of the observations, ~57” at 1.3 mm. The 
maps were made using uniform weighting. 
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lobe of the IRS3 CO outflow. The stellar mass of the IRS3 B source is ~8-12Mq, as reported by 
(2002). On the other hand, IRS3 A has a more extended IR nebulosity in the 


Preibisch et al. 


south-west direction. Although our CO data do not have sufficient spatial resolution to pinpoint 
whether IRS3 A or B is the driving source, the orientation of these nebulous features provides 
clues. Blue-shifted lobes pointing in the observer’s direction suffer less from extinction than the 
red-shifted lobes and the former are typically detected in the IR. Since the-direction of the IR 
nebulosity associated with IRS3 B (i.e. the microjet reported by Preibisch et al. 20021 ) matches 


the orientation of the blueshifted lobe observed with the SMA, IRS3 B is the more likely driving 
source of this outflow. 


Besides the IRS5 and IRS3 CO outflows, the large-scale redshifted outflow labeled D in the 
right panel of Figure |4] is closely aligned with an IR star detected in the NW corner of the image. 
This star is likely to be the outflow driver, consistent with the hypothesis that the blueshifted 
counterpart lies outside of the held of view to the northwest. The elongation of the knots on either 
side of IRS1 point to the location of the massive star IRS1 SW. Since this object is the exciting 
source of the MonR2 HII region (and it is thus at a late stage of evolution), it is possible that a 
low-mass companion is responsible for the CO molecular outflow labeled E. For the blue-shifted 
condensation close to IRS3 visible between -5 and -1.9 km s^ 1 (outflow F), IR observations reveal a 
nearby star aligned with the elongation axis, which could be its driving source. The exciting object 
of the blueshifted outflow C is less clear. The elongation axis suggests that one of the stars in the 
surroundings of IRS2 could be the driving source of outflow C, as indicated by the arrow in Figure 
[Tj For the other five additional outflow candidates (most of them located outside of the primary 
beam of the observations), their emission peaks appear associated with IR sources. The presence 
of IR stars at approximately the outflow geometric centers suggests that these outflows are true 
detections. 


3.2.3. Physical properties 


In this section, we estimate the physical properties of the six brightest CO outflows detected 
in MonR2 which present robust identifications (i.e. outflows A to F; see Figure [3J. 13 CO(2-l) 
emission was detected only in the blueshifted lobes of outflows A and B. The 12 CO emission from 
other outflow components without 13 CO detection was assumed to be optically thin. For the blue 
lobes of outflows A and B, where 13 CO(2-l) emission is detected, the optical depths are derived 
using Equation 1 in Roberts et ah. (201CJ) assuming a 13 C/ 12 C isotopologue ratio of ~l/69 (see 
Wilson 1999|). The derived values are 3.7 for the blue lobe of outflow A, and 11 for the blue lobe of 
outflow B. From this, the masses of these two lobes are higher by a factor of ~4 and 11, respectively. 
The correction factor is calculated as r/( 1 — exp [—r]), where r is the optical depth. 


We apply the method described in iGarden et al. (jl991f ) and IScoville et al. (jl986T ) to the 
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C0(2-l) transition: 


N C o = 1-08 x 10 13 (T ex + 0.92) exp ^ J T sdv cm 2 (2) 

Here 16.6 K is the upper level energy of the transition, dv is the velocity interval in km s' 1 , and 
T ex and Tb are the excitation temperature and brightness temperature, respectively. The outflow 
mass, momentum, energy, dynamical age, outflow rate and momentum rate are then given by: 


M = d 2 p g m(H 2 ) [ N CO (fi')dfi' 

(3) 

P = Mv 

(4) 

E = ^Mv 2 

(5) 

Length 

^dyn — 

^max 

(6) 

■ M 

M out = -— 

Myn 

(7) 

F t - P 
^dyn 

(8) 


(9) 


We assumed a CO to H 2 abundance of [CO/H2]=10~ 4 and a distance d of 830 pc (Herbst & Racine 


1976 ). p, g is the mean atomic weight of the gas (/x g =1.41) and m(H 2 ) is the mass of the hydrogen 


molecule. The momentum and energy associated with each velocity channel v in the ranges pre¬ 
sented in Figure [4] are added up to give total estimates for each outflow. The physical parameters 
of the CO outflows are reported in Table [5j Sin ce the shocked CO gas in outflows typic ally presents 
temperatures ranging from 30 K to 60 K (e.g. iBeuther et al. 1120021 : lArce et al. 1201(1 ) . we provide 
estimates of the physical parameters based on this temperature range. The masses we obtain are 
in the 10“ 4 to 10“ 2 M 0 range. The momenta are of order 10“ 3 - 10" 1 M 0 km s' 1 , while the energies 
are of order ~0.01-3.0 M 0 km 2 s" 2 . From these parameters, the dynamical ages of the outflows 
are found to lie between ~1000 and ~6000 yr, with corresponding outflow mass loss rates lying 
between ~0.1-17xl0" 6 M 0 yr" 1 . 


The out flow mass, momentum and energy are smaller than those from protostars of similar 


luminosities ([Zhang et al. 


2001 


2005|; IBeuther et al. 1 120021 ) . However, besides the optical depth 


correction, several additional effects could contribute to_ increasing our estimates: i) outflow gas 


at ve locities near the ambient clump velocity (lArce et al. 112010m ii) projection effects (lArce et al. 


2010): iii) 50% of the gas in an outflow could be in atomic form (R eipurth fe Ballv 200 1 j: iv) 


presence of high-velocity gas below the sensitivity of our observations ({Dunham et al. 2014); and 
v) the values presented here represent lower limits since the spatially extended emission in the 
outflow, which can be a large fraction of the total, may not be recovered in the SMA observations 
(e.g. Zhang et al., 2000; Qiu et al., 2009). 
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Table 5. Physical parameters of the CO outflows detected toward MonR2 with the SMA. 


Lobe 

Mass 
(10- 3 M 0 ) 

Momentum 
(10' 2 Mq km s" 1 ) 

Energy 
(M 0 km 2 s“ 2 ) 

^dyn 

(10 3 yr) 

Afout Pout 

(10- 6 M 0 yr 1 ) (10“ 5 Mq km s' 1 yr 1 ) 

A - IRS3 outflow 






Blue 

20.7-30.3 

25.2-38.0 

1.7-2.5 

1.8 

11.5-17.0 

14.1-20.7 

Red 

1.8-2.6 

2.8-4.2 

0.3-0.4 

2.0 

0.9-1.3 

1.4-2.1 

B - IRS5 outflow 






Blue 

4.4-5.5 

6.6-8.8 

0.55-0.77 

1.2 

3.6-4.6 

5.5-7.3 

Red 

1.6-2.4 

3.4-5.1 

0.4-0.6 

4.9 

0.3-0.5 

0.7-1.0 

C - Westernmost isolated outflow between IRS3 and IRS5 




Blue 

2.6-3.9 

3.1-4.6 

0.20-0.29 

1.5 

1.7-2.6 

2.1-3.1 

D - Large-scale redshifted lobe west of IRS5 




Red 

7.2-10.6 

8.3-12.2 

0.5-0.7 

6.2 

1.2-1.7 

1.3-2.0 

E - IRS1 outflow 






Blue 

1.6-2.7 

1.5-2.4 

0.07-0.1 

3.9 

0.4-0.7 

0.4-0.6 

Red 

0.2-0.3 

0.24-0.30 

0.01-0.02 

2.0 

0.11-0.13 

0.12-0.15 

F - Easternmost isolated outflow between IRS3 and IRS5 




Blue 

1.3-1.8 

1.2-1.8 

0.06-0.09 

2.1 

0.6-0.9 

0.6-0.8 


Note. — Estimates in columns 2, 3, 4, 5 and 6 are provided for a gas temperature range of 30 to 60 K. For 
the red lobe of B, the IRS5 outflow, the estimates include the knot up far to the northeast of IRS5 toward 
position (13”, 40”) in the left panel of Figure [5] (best visible in the primary-beam corrected map). Without this 
additional bullet, the mass, momentum and energy estimates are smaller by about 10%. 
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3.2.4- Spectral Energy Distribution (SED) fitting of IRS5 


To obtain further information about the physical properties of IRS5 and about its evolution¬ 
ary stage, we have built the SED of this source by using the fluxes measured by the SMA in 
compact configuration (Table Ep) and those observed in the near- and mid-IR, available in the lit- 


erature f Hackwel 


Gutermuth et al. 


et al. 


2009; 


de Wit et al. 


1982 Carpenter et al. Il997l : iKraemer et al. l 200 ll : iMueller et al. 120021 : 


.2009)- These values have then been fitted by means of the 


radiative transfer models of Robitaille et al. (120061 ). Note that like most of the IR observations, 
our data from the SMA compact configuration do not resolve the IRS5 A and B sub-sources iden¬ 
tified in VEX configuration (see Figure [2j) and therefore the SED fitting has been performed for 
source IRAS5 A+B. These results are reported in Figure [6l The best fit to the IRS5 SED provides 
a stellar mass of ~7M 0 and a stellar lumi nos ity of ~300Lq. This stellar luminosity is consis¬ 
tent with that previously derived by Hackwell et al. (1982) from 10-20 micron colors. The source 
presents a stellar age of ~ 10 1 yr, which is consistent with the idea of a young stellar object that is 
still accreting gas. We however point out that thi s es timat e is only approximate since the stellar 
evolutionary tracks used in the models of [R obitaille et al. (1200611 are fro m canonical non-ac c reting 

^1 996l: 


pre-main-sequence stars (Rernasconi & Maeder 


Siess et al. Il200fll h iRobitaille et al. (120061 ) 


indeed prefer to use an evolutionary classification into Stage 0/1, II and III sources instead of 
using stellar ages directly. Our source_ would therefore fall in Stage 0/1 after comparison to the 
SED shapes shown in Figure 7 of Robitaille et al. (2006). The IRS5 source is deeply embedded 
in the cloud with derived circumstellar extinction of A v ~50 mag, roughly consistent with the 
extinction inferred from the H 2 column densities measured toward IRS5 (i.e. ~20-40 magnitudes; 
see Section r3.3.2l) . The estimated disk mass is high (~O.2M 0 ), supporting the idea that IRS5 is 
relatively young. The derived disk accretion rate is ~10 ~ 6 M 0 yr' 1 , which is similar to the outflow 
rate calculated for outflow B (see Table [5]). We note that the disk parameters derived from the 
SED fitting should be taken with caution, due to the multiple nature of the IRS5 source and the 
relatively low-angular resolution (arcsecond scales) of the IR observations. The rich molecular line 
emission detected toward IRS5 (see Section 13.31) also indicates that IRS5 is at an early stage of 
evolution characterized by a hot-core-like chemistry. 


3.3. Molecular line emission from Other Molecules 

3.3.1. Morphology of the Molecular Gas toward MonR2 

In addition to the continuum emission 2 mm and 0.85 mm and to the 12 CO J= 2—M and 13 CO 
<7=2—M rotational lines, our SMA data report the detection of several molecular line transitions 
towards the MonR2 star forming region. The detected rotational lines are shown in Table O and 
include lines from molecular species typical of hot core chemistry such as CH 3 CN, CH 3 OH and 
sulphur-bearing molecules such as SO 2 , CS and OCS. The integrated intensity images of some 
representative lines are presented in Figure [3 
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Fig. 6.— SED fitting of the IRS5 source. Filled circles show the fluxes measured with the 
SMA (T able [2]) and those observed in the near- and mid-IR. which are a v ailable in the litera- 


ture (see 

Hackwell et al. 

1982; 

Carpenter et al. 

1997; 

Kraemer et al. 

2001; 

Mueller et al. 

2002; 

Gutermuth et al. 

2009; 

ie Wit et al. 

2009). Filled triangle indicates an upper limit to the 

mea- 


sured flux. The black line shows the best fit to the data while the grey lines indicate the 100 best 
fits to the IRS5 SED. The dashed line shows the emission arising from the star’s photosphere. The 
results of the SED fitting provide a stellar mass and luminosity of ~7 Mq and ~300Lq respectively 
for the IRS5 source. 
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Fig. 7.— Images of the MonR2 region for different spectral lines, indicated in the upper left-hand 
corner of each panel. The maps show the integrated emission from 5 to 15 km s' 1 , with each line’s 
central velocity set at 10 km s' 1 . As before, crosses mark the positions of the main IRS sources 
in the region and dotted contours indicate the negative 3 a level. The first contour and step level 
are 3a, except for OCS(19 -> 18), S0 2 (11 -> 10), CH 3 CN (K = 0 and K = 2), HCN(4 -»• 3), 
CS(7 —>• 6) and SO(8 —> 7), where the contour spacing is 6a. The beam size is shown as a red 
ellipse in the lower left corner: 2.99" x 2.98", P.A.= 86.6° at 1.3 mm, and 2.44" x 1.57", P.A.= 
114.7° at 0.85 mm. The dashed circle indicates the primary beam of the observations, ~ 57" at 
1.3 mm, and ~ 35" at 0.85 mm. The maps were made using natural weighting. 
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From Figure [TJ we find that the molecular emission in the MonR2 region i s concent rated to- 
ward the IRS3 and IRS5 sou rces. IRS3, with a luminosity of 1.5xl0 4 Lq (IHerming et al. Ill992l : 


Giannakopoulou et al. 1997), is expected to show a rich chemistry as a consequence of the evapo¬ 


ration of the mantles of dust grains as revealed by the detection of warm HtO, SO 2 and CH 3 OH 
molecular gas (excitation temperatures of 110-125 K for SO? and CH 3 OH and of 25 0 K for H 2 O; see 


van der Tak et al. 1 120031 : iBoonman fc van Dishoeck 1120031 : iBoonman et al. 1 120031 ) . However, only 


lines from SO 2 and the K= 0 to A'=5 transitions of CH 3 CN are measured toward this source. This 
is probably due to the fact that IRS3 falls outside of the primary beam of our observations, pre¬ 
venting the detection of fainter molecular lines. New high angular resolution observations centered 
on IRS3 are needed to study in detail the chemistry of the molecular gas toward this source. 


Figure [7] also shows that most of the molecular emission in our SMA images arises from IRS5 
and thus, besides IRS3, this source appears as the most chemically active object in the region. 
This confirms IRS5’s nature as a very young object with a chemically rich envelope. The molecular 
line emission toward IRS5 is very compact and its line profiles show a single Gaussian component 
with line widths ranging from 5 to 7 km s' 1 . We do not detect any broad velocity emission above 
the 3cr level in our images that could be associated with the high-velocity CO outflow reported in 
Section T3.21 This suggests that the molecular line emission detected toward IRS5 likely arises from 
the envelope surrounding this source, instead of from outflowing shocked gas. 


We note that the other two molecular condensations seen toward the north and south of IRS5 
in HCN, CS and SO, are side lobes generated by the limited UV plane coverage of our 0.85 mm 
observations. The presence of significant sidelobes in these images is due to a dynamic range 
problem as a result of the bright emission of these lines. 


3.3.2. Excitation Temperatures, Column Densities and Molecular Abundances toward MonR2 

IRS5 and IRS3 


Methyl cyanide 


Within our dataset toward MonR2, we have detected seven rotational transitions of the methyl 
cyanide molecule (CH 3 CN.) This molecule is commonly used to derive the excitation temperature 
of warm/hot gas in star-forming regions. In Figure 0 we report the images for the CH 3 CN K = 0, 
K = 2 and K = 4 transitions detected with the SMA. The CH 3 CN line spectra at the position 
of the IRS5 continuum peak were extracted from the image cube after primary-beam correction. 
The line parameters listed in Table [ 6 ] were derived from the fitting of the CH 3 CN K = 0 to K = 6 
transitions with single component Gaussian line profiles. The K = 0 and K = 1 transitions are 
very close in frequency and hence the accuracy in the fits suffers from partial merging of the lines. 
The exci tation temperature of the CH 3 CN gas was the n estimated using the rotational diagram 
method (jZhang et al. 1119981 : iGoldsmith fe Langer Ill999l . see also Appendix A). 


























Table 6. Derived parameters of the molecular line emission measured toward IRS5. 


Line 


v (GHz) 

E U [K ) 

Line Flux (Jy/beam km s J ) 

vlsr (km s x ) 

Av (km s x ) 

Ipeak (Jy/beam) 

CHgCN (12 -s- 11, K = 

0 ) 

220.74726 

69 

6.21 ±0.05 

9.34 ±0.08 

7.3 ±0.2 

0.801 ±0.018 

CH 3 CN (12 -» 11, K = 

1 ) 

220.74301 

76 

3.22 ±0.05 

9.4 ±0.3 

6.1 ±0.4 

0.492 ± 0.021 

CH 3 CN (12 -s- 11, K = 

2 ) 

220.73026 

98 

4.28 ± 0.06 

8.89 ±0.18 

6.9 ±0.4 

0.583 ± 0.020 

CH 3 CN (12 -s- 11, I< = 

3) 

220.70902 

133 

4.85 ± 0.06 

8.75 ±0.14 

7.4 ±0.3 

0.619 ±0.021 

CH 3 CN (12 -s- 11, K = 

4) 

220.67929 

183 

1.87 ±0.05 

8.89 ±0.21 

6.5 ±0.6 

0.268 ± 0.020 

CH 3 CN (12 -s- 11, K = 

5) 

220.64108 

248 

1.18 ±0.04 

9.20 ±0.15 

4.5 ±0.4 

0.247 ± 0.020 

CH 3 CN (12 -s- 11, K = 

6 ) 

220.59442 

326 

1.31 ± 0.06 

9.00 ± 0.30 

6.5 ±0.7 

0.190 ±0.022 

CH 3 CN (12 -s- 11, K = 

7) 

220.53932 

419 

< 0.14 

- 

- 

< 0.05 

SO 2 ( 111 , 11 —DOo.io) 


221.96521 

60 

4.47 ±0.05 

9.09 ±0.14 

6.1 ±0.3 

0.646 ± 0.020 

OCS (19 -s- 18) 


231.06098 

111 

5.35 ± 0.07 

8.95 ±0.09 

6.5 ±0.2 

0.78 ±0.03 

13 CS (5 -s- 4) 


231.22099 

33 

3.14 ±0.15 

9.22 ±0.23 

6.2 ±0.5 

0.48 ± 0.06 

CH 3 OH (10 2 ,9 ^9 3 ,6 A 

“) 

231.28111 

165 

4.07 ±0.07 

8.76 ±0.11 

5.8 ±0.2 

0.66 ±0.03 

CH 3 OH (13i,i2 -> 13 0 ,i3 A) 

342.72980 

228 

20.6 ±0.3 

9.07 ±0.12 

6.1 ±0.3 

3.16 ±0.16 

CS (7 -1 6 ) 


342.88300 

66 

42.8 ± 0.23 

8.59 ±0.01 

3.72 ± 0.03 

10.8 ± 0.1 

HC 1S N (4 -s- 3) 


344.20032 

41 

6.8 ±0.3 

8.21 ±0.24 

5.2 ±0.6 

1.24 ±0.16 

SO (8 -1 7) 


344.31061 

88 





Component 1 




15.8 ±0.2 

7.3 ±0.7 

6.4 ±0.7 

2.3 ±0.1 

Component 2 




3.42 ± 0.18 

9.7 ±0.7 

4.4 ±0.7 

0.7 ±0.1 

HCN (4 -s- 3) 


354.50548 

43 

60.2 ±0.7 

8.18 ±0.03 

5.4 ±0.1 

10.4 ±0.4 


Note. — The errors on the line flux were calculated as ay/Sv x Av, where 5v is the velocity resolution in the spectrum and Av is 
the linewidth (column 6 in this Table). Since the CH 3 CN ( K= 7) transition was not detected, the upper limit on its peak intensity 
is given by the 3cr level in the spectrum, and the upper limit on the integrated intensity as 3 ay/5v x Av. The SO (8 —> 7) transition 
presented with a double Gaussian profile corresponding to two velocity components. 
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Fig. 8.— Methyl cyanide rotational diagram for the K = 0 to K = 6 ladder detected toward IRS5. 
The data are shown as black plus signs with corresponding lcr error bars, and the linear fit as a 
black line. The upper limit derived for the K = 7 transition is shown as a black arrow. 
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The rotational diagram of the CH 3 CN molecule, derived by following the method described 
in the Appendix, is shown in Figure [HJ From the slope of the linear least squares fit to the 
data, we estimate an excitation temperature for the CH 3 CN gas of 144 ± 15 K. This value is 
cons istent with typical gas temperatures found in massive ho t core sources such as Cepheus A HW2 
(see iMartfn-Pmtado et al. 1120051 : Ijimenez-Serra et al. 1120091 ). The derived total column density of 
CH 3 CN toward IRS5 is (3.2 ± 0.4) x 10 13 cm" 2 . 

Although IRS3 falls outside the primary beam of the SMA observations, we can provide an 
estimate for the excitation temperature of the gas toward this source from the measured CH 3 CN 
line emission. The deri ved e xcitat ion temperature is 126 ± 22 K.This temperature is consistent 
with those measured by van der Tak et al. (120031 ) from single-dish observations of several SOo and 
CH 3 OH lines (~110 K to 125 K). We note that we provide no additional details for the CH 3 CN 
rotational diagram calculated for IRS3 because the K > 2 lines are detected only at the ~ 3a level 
(the source falls at the edge of the primary beam). 

Other molecular species 

The column densities of the other molecular species detected toward IRS5 are estimated by 
assuming LTE and optically thin emission. In the case of HCN, the emission is not optically thin 
because HC 15 N is detected (see Table[ 6 j). Therefore, we correct for optical depth (see th e expr ession 
for the correction factor in Section l3.2.3l) based on an isotopic ratio of 1/388 (Wilson 19991)- The 
derived optical depth for HCN is r~47. We use Equation IA3I to determine the values of N u /g u 
for every rotational transition, and then calculate the total column density, Wot of the molecular 
species by reformulating lA2l into: 


Wot = —Q(T mt )e E ^ kT ™K 
9u 


( 10 ) 


The values of Sp? and Q(T m t — 150 K) are taken from the Cologne Database for Molecular 
Spectroscopy (CDMS; Mueller et al. 112005 b We estimate the molecular abundances, %, by matching 
the beams of the 1.3 mm and 0.85 mm continuum maps with those of the line maps at 1.3 mm 
and 0.85 mm, respectively. The H 2 column density toward IRS5 is calculated from the measured 
1.3 mm and 0.85 mm peak fluxes of IRS5 (Table [2]) following the method in [Enoch et al. 1 ( 2006 ). 
We assume a gas-to-dust ratio of 100. The derived H 2 column densities toward this source are 
1.7 x 10 22 cm ' 2 from 1.3 mm and 3.8 x 10 22 cm ' 2 from 0.85 mm. These values are derived from 
the peak fluxes at different angular resolutions to match the peak flux of the molecular lines. The 
resulting molecular column densities and abundances are reported in Table [71 























Table 7. Molecular column densities and abundances measured toward IRS5 with the SMA. 



ch 3 cn 

S0 2 

OCS 

13 CS 

CH3OH 

CS 

hc 15 n 

SO 

HCN 

N (x 10 15 cm" 2 ) 
X (xHT 8 ) 

0.032 ± 0.004 

0.2 

1.23 ±0.11 

7.3 

1.31 ±0.08 

7.7 

0.031 ± 0.003 

0.2 

24.5 ±2.0 

150 

46.8 ± 0.004 

1.2 

0.027 ±0.001 

0.07 

1.12 ±0.02 

3.0 

9.87 ±0.001 

26.3 


to 
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Discussion 


4.1. Outflow feedback in the MonR2 star forming cluster 


The presence of multiple outflows from a cluster of low- to intermediate-mass stars may have 
important implications in cluster and star forma tion, since outflows can inject large amounts 
of turbulence into the surrounding environment (|Li &; Nakamura 120061 : iNakamura fc Li 120071 : 
Carroll et ah 2009)- In this section we discuss the impact of the newly discovered cluster of com¬ 
pact CO outflows in the center of the MonR2 clump, comparing with the one due to the already 


known large-scale outflow (IWolf et al. 1 11990). 


The outflow feedback in a given environment depends on the outflow momentum as well as 
the mass a nd radius of the region being affected. The outflow momenta are presented in Tabled 
Ri dge et ah (120031 1 estimated a mass of ~1800 M 0 for the MonR2 clump within a radius of ~0.7 


pc 


Walker et ah 


fl 1990 1 found that the inner ~0.22 pc of the clump contains ~180 Mm. These 
values are consistent with the density profile of r^ 1 found in the region by Choi et ah (2000). 


The cluster of outflows presented in this paper is distributed within the inner ~0.15 pc radius of 
the clump. Scaling the mass using the r _1 density profile resulting from the estimates above, we 
obtain a mass of ~75 M 0 for this region. In the following analysis, we will use R reg = 0.15 pc and 
M reg = 75 M 0 to evaluate the impact of the cluster of com pact S MA outflows. In the case of the 
CO giant outflow, whose size is much larger (~ 0.6—7 pc; Wolf et ah 19901 ) — and henc e acts o ver 


a Jarger and more massive region— we will use the values of the clump studied by IRidge et ah 
(H): Rreg = 0.7 pc and M reg = 1800 M 0 . In Table [8] we summarize the physical properties 


xeg — 1 ctixvj. j. vj . re g 

of the two regions of interest, along with the outflow feedback parameters to be explained in the 
following subsections. 


4-1.1. Turbulent support of the MonR2 clump 


In the absence of replenishment of the initial turbulence, the MonR2 clump is expected to 
collapse in a free-fall time, t& The free-fall time for the region with R ieg = 0.15 pc is estimated 
to be iff ~lxl0 5 yr, and for R rC g= 0.7 pc it is iff ^2xl0 5 yr. However, the estimated average age 
of the stellar cluster is larger than these values, ~1 Myr (Carpenter et al. 19971; lAndersen et ah 


2006), which implies that a mechanism may have pre vented the collapse of the MonR2 clump in 


its free-fall time, allowing the formation of stars. Wolf et ah (jl990i ) already claimed that the 


laps 

3 ) 


large-scale CO outflow detected in the region could provide enough support to the whole clump. 
We evaluate here whether the population of CO outflows revealed by the SMA can also play a 
significant role. 


By following the analysis of Rivilla et ah (2013b), we compare the rate at which the initial 


O / O 3M 

ts = 32 q p , where G is the gravitational constant and p is the gas density, calculated as p = 47rfl r a et 
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turbulence decays, L tur b, with the rate at which the clump gains energy due to outflows, L ga j n . 
The former is given by L tur b = RturbAff) with -E tul -b = (3/2)M reg af D , where ctid is the observed 
one-dimensional velocity dispersion of the gas in the clump and M rPP the_clump mass. We evaluate 


the ID velocity dispersion using the CS linewidth observed by Tafalla et ah_ (1997]). These authors 


showed that the linewidth decreases with the radius. Prom their Fig. 14 we consider an average 
value of linewidth of Av ~2.3 km s^ 1 for radius < 0.15 pc and ~1.9 km s -1 for radius < 0.7 pc. 
For a ID Maxwellian distribution of velocities, the relation between the linewidth and the one¬ 
dimensional velocity dispersion is Av = 2\/2In2 ctid, and hence we obtain ctid ~0.96 km s -1 and 
0.81 km s -1 , respectively (Tabled]). 


On the other hand, L ga i n is given by 


£ core ^3D 


= ^^outO'iD, where M C( 


core is the mass set 

in motion by outflows per unit time and F out =P ou t/tdyn is the outflow momentum rate. From 
the measured values of the CO molecular outflows (see Tabled]) we obtain that T ga m/Tturb ~0.3 
(Table[8]). A s sh own in Figured] this value is similar to the ones found in several Perseus clumps 


(Arce et al. 2010). 


We note that the values of the outflow para m eters d erived from the SMA data have likely been 
underestimated because of several reasons. Arce et al. (2010j) suggested correction factors for the 
outflow momentum of 1.4 and 2 to account for inclination effects and for outflow gas at velocities 
near the ambient clump velocity. In Figure [9] we have included a ” correction vector” to account 
for these two effects, which yield a T ga ; n /L tur b ratio of ~ 0 . 9 . This value is a lower limit, because 
the gas in an outflow could be in atomic rather than molecular form due to dissociative shocks 
(jReipurth &; Bally 2001), and because some high-velocity emission from the outflow may be below 
the sensitivity of our observations. 


Recently, Dunham et ah 2014 quantified these corrections in a sample of single-dish outflow 
observations. Here we have used the mean corrections factors found in this study (factor of 41 
and 93 for ^outflow and ^outflow) respectively). Moreover, we have considered an additional fac¬ 
tor of 2 to account for the presence of missing extended flux due to short-spacings not covered 
by the interferometer. Figure [9] therefore also includes a second ’’correction vector” to account 


for all these effects. The uncorrected data from outflow observations by lArce et al. (J2010|) and 
St anke fe William s (2007]) are shown for comparison. Considering all the corrections, the ratio 
A ga in/Aturb for the SMA outflows would reach a value of ~61 (Fig. [9]). This implies that even 
if only a fraction of the energy in the compact outfows is transformed into turbulent energy, the 
SMA cluster of outflows could help to maintain the turbulence in the MonR2 clump center. For 
the large-scale CO outflow, we obtain L ga m/Dturb ~122 after applying all corrections except for the 
missing flux correction. 
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4-1.2. Broadening of the linewidths induced by the population of the CO outflows 


As we already mentioned, Tafalla et aL ([19971) found that the linewidths of the dense molecular 
tracer CS increase towards the center of the MonR2 region. They suggested that this is due to 
a combination of the broadening produced by the giant CO outflow and an additional turbulence 
source in the inner region. Since we have detected a cluster of young molecular CO outflows at 
the center of the MonR2 clump, we evaluate whether the outflow feedback can contribute to the 
broadening of the molecular linewidths. 

Assuming that the total momentum of the SMA outflows is transferred to the gas, we obtain 
that the ratio between the velocity dispersion produced by outflows cr ga i n = ° w and the observed 
velocity dispersion ctid is av a . m /o m « 1 ( see Table [ 8 ] and Figure [9]). However, considering the 
correction factors found by Dunham et ah (2014), this value reaches <T ga i n /criD ~ 0.9. This implies 


that the turbulence injected by the cluster of CO outflows in the central region of MonR2 could be 
responsible for the observed velocity dispersion. In the case of the large-scale outflow, the corrected 
ratio Ugain/cJiD is >>1 (see Tableland Figure [9]), indicating that this agent may also contribute 
significantly to the velocity dispersion observed in the whole MonR2 clump. 

The turbulence in jected by th e com pact SMA outflows may also explain the difference in 
linewidth observed by Tafalla et al. (ll997n between the center and the periphery of the clump. This 
difference is Audifr ~0.5 km s -1 ( Tafalla et al. 19971 ). which is equivalent to a velocity dispersion of 
(Jdiff = 0.21 km s _1 . This value is lower than our corrected value for cr gain ~ 0.85 km s , suggesting 
that the broadening of the molecular linewidths in the center of the cluster may be produced by 
the compact CO outflows. 


4.2. IRS5: An intermediate- to high-mass star in the hot-core phase 


In Section lT3l we have reported the detection of rich molecular line spectra toward the IR 
sources IRS3 and IRS5 in the MonR2 cluster. These spectra present several molecular line transi¬ 
tions from complex organic species, such as CH 3 OH and CH 3 CN, and from sulfur-bearing species, 
such as SO 2 or OCS. The high temperatures measured toward the innermost envelope regions 
around massive stars (of few 100 K) cause the thermal evaporation of the icy mantles of dust 


hot cores ; see e.g. Charnley et al. 


1992 


Charnley 

1997; 

Yiti et al. 

2004; 

Wakelam et al. 

2004) 


For IRS3, previous millimeter and infrared observat i ons toward this source a lready noted a rich 
hot-core chemistry (jBoonman fc van Dishoeck 120031 : Ivan der Tak et al. 1200.31 ) . Although shocks 
are also expected to enhance the abundance of certain molecul ar spe cies in star forming regions, the 
radiative transfer modelling study carried out bv lvan der Tak et al. (1200.31) showed that their con¬ 


tribution is likely very little. IRS3 is a cluster of IR sources (jPr eibisch et al . 200 2) and therefore 
the most massive objects in the region, IRS3 A and B (with derived star masses of ~1O-15M 0 ), are 
the best candidates for driving the hot-core chemistry found toward the IRS3 cluster. Higher angu- 
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Fig. 9.— Outflow feedback from the MonR2 cluster of outflows. The blue filled square corresponds 
to the SMA outflows presented in this paper (correction factors not in clude d). We compare this 


feedback with the impact of the giant outflow (magenta open star, Wolf et ah 1990|) and several 


clusters of outflows from other star-forming regions. The short black arrows in dicat e the correction 
accounting for the underestimation of the ou tflow parameters follow ing lArce et al. (|2010l ) . Longer 
arrows account for the corrections found by Dunham et al. (2 014 ) and an additional factor of 2 
due to the lack of short-spacings. 


lar resolution interferometric observations are needed to determine the origin of the rich molecular 
line emission arising from this source. 

For IRS5, our SMA observations reveal that this object also presents a hot-core-like chemistry, 
indicative of its youth. The molecular abundances measured toward IRS5 indeed range from 2x 10 ~ 9 
for CH 3 CN to 2xl0 ~ 6 for CH 3 OH, showing that they have been enhanced by several orders of 
magnitude with respect to the quiescent gas in molecular dark clouds, as expected in hot cores 
(see e.g. van der Tak et ah 200(1 200a) ■ Supporting this idea, the derived excitation temperature 
of the gas toward IRS5 is typical of hot cores (~150K, see Section l3.3.2lh indicating that the rich 
chemistry observed toward this source is a consequence of the evaporation of large amounts of 
molecular material from the icy mantles of dust grains. 

Sulfur-bearing molecules are considered good chemical clo cks si n ce t heir abunda nce a re pre- 
dicted to vary by several orders of magnitude with time (e.g. 


Charnlev 

1997; 

Viti et al. 

2004) 


Sulfur chemistry in hot cores starts with the injection of ELS into the gas phase after grain mantle 
evaporation (although H 2 S is undetected in solid form in the ISM, th is sp ecies is believed to be the 
most abundant carrier of atomic sulfur on the surface of dust grains; Smith 1991). H 2 S is destroyed 
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by H and H 30 + to form S and H 3 S + , which then react with O and OH to form SO (Charnley 1997). 
For temperatures <300 K, SO is destroyed in favour of SO 2 through the reaction O + SO —» SOo + 
photon, which t akes over f or time-scales on the order of 10 4 * yr, yielding a x(S 0 )/x(S 02 ) abundance 
ratio <1 (Wakelam et al. 2004). Keeping in mind the uncertainties in our abundance calculations, 
the abundance ratio x(SO)/x(SO2)~0.4 measured toward IRS5 implies the age of this source to be 
at least 10 4 yr. This chemical timescale is in reasonable agreement with the age estimated from the 
SED fitting (Section ^. 2. 21) . The H 2 densities are likely >10 6 cm~ 3 . si nce SO 2 is never found to be 
more abundant than SO for lower H 2 densities (Wakelam et ah 2004!) • The derived abundances of 
OCS and CS (~8xl0” s and ~10 -8 , respective ly) are also simi l ar to those predicted by chemical 
models for the same time-scales of >10 4 yr (see ICharnlev 1119971 : IWakelam et al. Il2004l l. 


4.3. Competitive accretion as possible origin of the MonR2 star cluster 


Previous near-IR observations carried out bv lHodapp (120071 ) toward the M onR2 star cluster re¬ 
ported a population of H 2 jets surrounding the MonR2 clump center. From this, Carpenter k Hoda pp 
(J2008) proposed that the observed distribution of H 2 2.12 ^urn outflows was consistent with a sce¬ 
nario where star formation had been triggered by the interaction of the large-scale CO outflow 
with the surrounding ambient cloud. Our SMA observations toward this cluster, although limited 
to the inner l'-region toward MonR2 (or ~0.25pc at 830pc), reveal a population of 11 new CO 
outflows that appear widely distributed across the central region of the clump. This population of 
CO outflows, which was not d etected in near-IR or mid-IR surveys mainly du e to the high extinc¬ 
tion toward the cluster center ( Carpenter et al. 1997 : Gutermuth et al. 2005 ^1. is likely associated 
with the youngest population of Class 0 sources in the cluster. Since this population is widespread 
across the MonR2 clump, an alternative mechanism different from triggering must dominate the 
star formation process in the MonR2 cluster. 

By inspecting the distribution of Class I and Class II sources in MonR2, Gutermuth (2005) 
found that the Class I sources cluster around the MonR2 region along the north-south direction 
while the more evolve d Class I I objects ar e distributed across the cluster (see Figure 13 in the 
review paper of lCarpenter k Hodapp 1120081 ) . In additi on, the distributio n of molecular gas, which 
was de rived from the multi-tra nsition analysis of CS (jChoi et al. 120001 ) and the column density 
of CO (Xie k Goldsmith 1994), peaks toward the center of the cluster. This configuration would 
be consistent with the competitive accretion scenario where the youngest (and the most deeply 
embedded) YSOs, probed by the CO outflows detected with the SMA, form at the center of the 
MonR2 cluster as a consequence of gas accretion by the potential well. This would explain not 
only the differences in the distribution of Class 0, I and II sources toward MonR2, but also the 
centrally-peaked gas distribution toward the cluster center. 


4 The Spitzer data also suffer from saturation and crowding toward the MonR2 cloud center (see [Gutermuth et al. 


200 - 3 ) - 
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If competitive accretion were the dominant mechanism for star formation in the MonR2 cluster, 
one would expect the most massive objects in the cluster (IRS 1, IRS 2, IRS 3 and IRS 5) to grow 
via gas accretion accompanied by a significant population of low mass stars (jBonnell & Bate 2006). 
As seen in Figure [4] the density of low-mass stars in MonR2 is found to be higher around the most 
massive objects in the cluster such as IRS1, IRS2 and IRS3. A detailed study of the distribution 
of the stellar density in MonR2 (Rivilla et al. 2015, submitted) indeed shows that this distribution 
is centrally peaked coinciding with the location of the IRS1, IRS2, IRS3 and IRS5 sources. The 
maximum value of the stellar density distribution (>10 5 5 stars pc -3 ; Rivilla et al. 2015, submitted) 
is found close to the location of the IRS2 source, which corresponds approximately to the geometrical 


comparable masses 


(Preibisch et al. 

2002; 

1 --— 

Jimenez-Serra et al. 

2013; 

Rivilla et al. 


2014). The presence of a few scattered sub-clusters hosting massive stars is an expe cted outcome 
of the competitive accretion scenario, as shown by the simulations of Bon nell et al. (2011). 


The competitive accretion scenario could also explain the large-scale CO outflow detected in 
Mon R2, which emanates from the cluster but for wh ich no powering source has been identified 


yet (ITafalla et al. [1997]; iGiannakopoulou et al. 119971 ) . Since the stellar densities at the center of 


the MonR2 clump are very high, close dynamical encounters are likely to occur as proposed for 


Zapata et al. 

2009, 

2013; 

Rivilla et al. 

2013a, 

2014) 


2014). A large-scale SMA 


mosaic covering the extent of the MonR2 CO outflow is needed to verify this hypothesis. 


5. Conclusions 

Our SMA study of the MonR2 clump center has characterized the thermal continuum and 
molecular line emission from the main sources in the region at (sub-)millimeter wavelengths. The 
main conclusions of this study are: 

• Continuum observations: IRS2 and IRS5 are the brightest continuum sources detected at 
1.3 mm and 0.85 mm towards MonR2 (IRS3 falls outside the primary beam of the observa¬ 
tions, preventing us from accurately characterizing its physical properties). While IRS1 and 
IRS2 are dominated by free-free emission (oirsi ~ —0.5 and «irs 2 ~ —0.2), IRS5 is domi¬ 
nated by dust thermal emission («irs 5 ~ 3.0.) We obtain a mass estimate of rs -' 0.1 — 0.3 M 0 
for IRS5, and a lower limit of ~ 0.05 — 0.15 M 0 for IRS3, the most luminous source in the 
region. Higher-angular resolution observations in the SMA Very Extended configuration re¬ 
veal the substructure of IRS5 and IRS1 at the 0.5” level. IRS2 however retains its compact 
morphology at this resolution. 

• Molecular outflows: We detect 11 previously unknown young outflows in the 12 CO (2-1) tran¬ 
sition. IRS5 and IRS3 are associated with bright collimated bipolar outflows, with masses 
in the 10' 3 -10" 2 M 0 range. Their momenta are of order 10" 2 - 10' 1 M 0 km s' 1 , while the 
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energies are of order ~0.3-2.5 Mq km 2 s“ 2 . A comparison with known IR sources indicates 
that the IRS5 outflow is likely associated with the IRS5-A source and that the IRS3 out¬ 
flow is likely drive n by the near IR source IRS3B, whose mass was reported to be 8-12 Mq 


(Preibisch et al. 2002). The SED fitting for IRS5 yields a stellar mass and luminosity of 


~ 7 Mq and ~ 300 L 0 , respectively, and further suggests that IRS5 is a deeply embedded 
young stellar object still accreting gas. 


• IRS5 chemistry: Our observations show that, besides IRS3, IRS5 is the most chemically active 
source in the region. From the methyl cyanide K ladder transitions, we derive excitation 
temperatures of 144 ± 15 K for IRS5 and of 126 ± 22 K for IRS3. We also provide estimates 
for the molecular column densities and abundances for the molecular species detected towards 
IRS5. 


We conclude that IRS5 is an intermediate- to high-mass star at an early stage of evolution, as 
shown by its rich hot-core-like chemistry and the bright, collimated molecular outflow originating 
from it. The CO outflows are found to contribute to the maintenance of turbulence in the MonR2 
clump, helping to support it against a global rapid collapse. Together, the newly discovered outflows 
might have significantly contributed to the velocity dispersion in the MonR2 clump and therefore 
in the broadening of the molecular linewidths toward the clump center. Finally, the detection of 11 
new outflows widespread across the MonR2 clump would be consistent with a scenario where star 
formation in the MonR2 clump takes place through competitive accretion as opposed to triggering. 
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number PIIF-GA-2011-301538. V.M.R. acknowledges financial support from the Spanish project 
AYA2010-21697-C05-01. 
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A. Rotational Diagram Method for the CH 3 CN rotational lines 


We detect seven rotational transitions of the methyl cyanide (CH 3 CN) molecule, commonly 
used to derive the kinetic temperature in star-forming regions, toward MonR2. Images for the 
CH 3 CN k = 0, k = 2 and k = 4 transitions are presented in Figure [71 In Table [ 6 ] we present the 
parameters extracted from fitting each K = 0 to K = 6 transition to a Gaussian emission line, as 
well as upper limits on the integrated and peak intensities for the undetected K = 7 transition. 
The K = 0 (220.74726 GHz) and K = 1 (220.74301 GHz) transitions are very close in frequency 
and hence the fits suffer from partial merging of the lines. 


We apply the rotational diagram method to CH 3 CN in order to estimate the gas excitation 
temperature of IRS5. To this end we assume the emission is optically thin, and that the medium is 
in Local Thermodynamic Equilibrium (LTE). This implies that the level populations are described 
by the Maxwell-Boltzmann distribution with temperature T ro t, such that: 


_ Su c —hv/kT I0 t (Al) 

Ni gi 

where g u and gi designate the statistical weights of the upper and lower levels respectively, and v 
is the transition frequency. Rewriting the expression after assuming LTE gives: 


Nu 

9u 


N tot ~E u /kT Ioi 

QCTrot) 


(A2) 


where Q(T) is the partition function, E u is the upper level energy, and N tQ t is the total column 


density. The quantity is calculated following Zhang et ah (1998): 


Njk 

9jk 


= 1.669 x 10 


17 f Tftdv(K km s' 1 ) 
S.j /c /i( D) 2 // (MHz) 


an 


-2 


(A3) 


where g is the dipole moment and equals 3.92 Debye (Blake et al. 19871); Sjk is the line strength 
and is given by Sjk = ( J 2 — K 2 ) / J (J = 12, K = 0 to K = 6 here); and Tb is the brightness 
temperature. The velocity-integrated brightness temperature can be obtained from the line areas 
though the Rayleigh-Jeans law: 


T B dv = 1.224 x 10 6 


S v 


zz(GHz ) 2 M 


"\2 


K km s 


-1 


(A4) 


where S v is the line area in Jy/beam km s " 1 as in Tabled and 0b is the geometric mean of the major 
and minor axis of the beam in arcseconds. The statistical weights are given by gj k = (2 J + 1 )gx 
where gx = 4 for K = 0 or I( / 3n, and gx = 8 for K = 3n, where n is a nonzero integer. 


The resulting rotational diagram is shown in Figure 0 Extracting the slope from a linear fit to 
the data following Equation IA2I gives an estimate of the gas excitation temperature of 144 ± 15 K. 
The error of T ex has been derived by p ropaga ting errors. This value is consistent with typical 
temperatures of hot cores (e.g. Jiinenez-Serra et al. 2009). The total column density can easily be 
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estimated from the y-intercept of the fit. The CH 3 CN partition function is given by Blake et al. 
d 19871 ): 


Q(Trot) 



(A5) 


A and B are the rotational constants, with A = 158099 MHz and B = 9198.9 MHz from the Cologne 
Database for Molecular Spectroscopy (CDMS; Mueller et al. 2005). This yields an estimate for the 
column density of (3.2 ± 0.4) x 10 13 cm" 2 . 
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Table 8: Sizes and masses of the regions where the SMA cluster of compact outflows and the 
large-scale outflow are considered to act, along with the derived outflow feedback parameters. 


Outflow 

R Y Qg 

(pc) 

M reg 

(M 0 ) 

T/turb/-hg ain 

^ID/^gain 

SMA cluster 

0.15 

75 

0.3 (0.9) [61] 

0.01 (0.03) [0.9] 

large-scale CO 

0.7 

1800 

0.7 (1.8) [122] 

0.7 (2) [61] 


The values in parentheses a re corrected for gas at ambient velocities and projection effects 


following Arce et ah (2010). The values in square brackets are corrected for gas at ambi¬ 


ent velocities, projection eff ects, atomic/molecular gas fraction, sensitivity (using the mean 


correction factors from Dunham et ah 2014). For the cluster of outflows observed with the 


SMA, the effect of missing flux was also corrected (see Section 14.1.11) . 
















